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ABSTRACT 


Fluid  production  from  both  normally  and  abnormally  pressured 
reservoirs  has  long  been  recognized  as  a  function  of  the  formation 
fluid  properties  and  the  reservoir  rock  parameters.  Fluid  production 
itself,  however,  changes  the  stress  state  of  the  reservoir  resulting 
in  rock  deformation  which  affects  total  fluid  and,  therefore,  energy 
recovery. 


This  study  was  undertaken  as  part  of  a  project  to  determine 
and  predict  the  effects  of  fluid  production  on  the  Pleasant  Bayou  geo- 
pressured,  geothermal  reservoir  in  Brazoria  County,  Texas.  Data  was 
collected  from  triaxial  loading  tests  conducted  on  sandstone  samples 
cored  from  the  potential  producing  horizons  of  the  reservoir.  All  tests 
were  conducted  at  overburden  and  pore  pressures  which  simulated  insitu 
conditions.  Empirical  relations  were  developed  which  reflected  changes 
in  the  resistivity,  porosity,  and  permeability  as  functions  of  the 
effective  stress.  Relationships  among  the  parameters  were  found  to  be 
of  practical  use  in  predicting  formation  changes  as  the  reservoir  stress 
state  changed  due  to  fluid  withdrawal. 

Results  showed  that  the  resistivity,  porosity,  and  permeability 
were  dependent  on  the  effective  stresses  in  the  reservoir.  The 
cementation  factor  in  Archie's  equation  was  also  found  to  be  a  pressure 
function.^ In  contrast,  using  the  electrical  and  hydraulic  conductivity 
analog,  the  product  of  the  formation  factor  and  the  permeability  was 


found  to  be  Independent  of  pressure. 


Page 
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I.  INTRODUCTION  AND  SCOPE  OF  THE  INVESTIGATION 


1.1  Introduction 


With  alternate  forms  of  energy  becoming  more  and  more  economically 
attractive,  much  interest  and  research  is  b<Hng  focused  on  geopressured , 
geothermal  reservoirs.  These  reservoirs  contain  significant  amounts  of 
methane  held  in  solution  and  are  air*  a  viable  thermal  energy  source. 
However,  fluid  production  over  time  -  i.  .urt  jes  the  stress  state  of  the  for¬ 
mations  which  car.  detrimentally  affect  the  ultimate  energy  recovery. 
Therefore,  a  broad  knowledge  of  the  deformation  characteristics  of  the 
reservoir  rocks  and  of  the  resulting  changes  in  permeability,  porosity, 
resistivity,  and  bulk  volume  is  a  desirable  goal. 

Changes  in  the  permeability,  porosity,  and  resistivity  which  re¬ 
sult  from  changes  in  the  stress  state  have  been  approached  from  several 
directions.  They  have  been  related  to  deformation  of  rock  under  pressure 
and,  therefore,  have  been  extensively  studied  as  direct  functions  of 
pressure.  Many  investigators  have  approached  these  parameters  from  this 
"rock  mechanics"  point  of  view.  Most  of  the  studies,  however,  have  been 
conducted  under  atmospheric  pore  pressure  conditions  which  do  not  fully 
simulate  the  lnsitu  state. 

Another  approach  has  been  to  relate  the  resitivity  of  rock  to  its 
permeability  and  porosity.  This  approach  is  quite  practical  since  the 
formation  resistivity  is  a  measurable  parameter  via  electric  logging. 

The  present  work  represents  both  of  these  viewpoints  and  was 


undertaken  to  study  the  Influence  of  the  effective  stress  on  the  re- 
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sistivity,  porosity,  and  permeability  of  samples  from  the  Gulf  Coast 
geothermal  province.  Specifically,  an  attempt  has  been  made  to  deter¬ 
mine  the  empirical  relationships  among  these  three  parameters  and  be¬ 
tween  them  and  the  effective  stress  at  elevated,  insltu  pore  pressures. 
Tests  were  conducted  on  samples  extracted  from  wells  drilled  in  the 
Department  of  Energy's  Pleasant  Bayou  Geothermal  Reservoir  in  Brazoria 
County,  Texas.  Figures  1  and  2  depict  the  province  and  reservoir 
locations . 

The  Pleasant  Bayou  prospect  is  characterized  by  several  geo- 
pressured,  geothermal  zones.  Pressures  exceed  the  normal  hydrostatic 
saltwater  head  by  150%  and  correspond  to  a  pressure  gradient  as  high  as 
0.78  psi/ft.  Temperatures  range  from  250°F  to  350°F.  According  to  Be- 
bout  et  al.1  the  reservoir  with  the  best  overall  quality  is  a  medium  to 
coarse-grained,  bedload,  fluvial,  channel  deposit  at  14,640  to  14,710 
feet.  These  Frio  reservoir  sands,  however,  are  also  found  in  other 
distinct  zones  from  11,750  to  16,050  feet.  Samples  were  extracted  from 
various  depths  within  this  depth  range. 

This  work  was  undertaken  to  obtain  experimental  data  to  predict 
the  effects  of  fluid  production  on  the  reservoir.  As  the  reservoir  is 
exploited,  formation  pressures  decline  resulting  in  an  increase  in  ef¬ 
fective  stress  and  causing  possible  formation  compaction,  porosity  and 
permeability  reduction,  and  a  fluid  recovery  decrease. 

Empirical  equations  were  developed  based  on  experimental  results 
and  previous  works1 1" 1 5* 1 7* 21  * 2S  which  showed  the  interrelationships 
among  the  porosity,  permeability,  and  resistivity.  These  equations  can 
be  practically  applied  and  are  useful  in  estimating  the  changes  in 
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porosity  and  permeability  with  effective  stress  throughout  the  producing 
life  of  the  reservoir. 


1.2  Scope  of  Investigation 

A  review  of  literature  shows  that  most  previous  studies  on  re¬ 
sistivity  measurements  have  been  conducted  under  atmospheric  pore  pres¬ 
sure  conditions,  which  do  not  adequately  simulate  the  insitu  stress 
state.  In  addition,  none  of  these  investigations  have  resulted  in  prac¬ 
tical  equations  for  use  in  estimating  changes  in  the  resistivity,  per¬ 
meability,  and  porosity  with  decline  of  formation  pressure. 

The  objective  of  this  work  was  to  formulate  empirical  relation¬ 
ships  for  these  parameters  as  functions  of  the  effective  stress.  For¬ 
mation,  overburden,  and  confining  pressures  were  applied  to  the  rock 
samples  which  essentially  duplicated  what  were  considered  to  be  insitu 
conditions.  Pore  pressures  were  reduced  to  simulate  fluid  production, 
and  measurements  were  taken  of  the  resultant  resistivity,  permeability, 
and  porosity  changes.  Results  were  used  to  check  the  validity  of  the 
Archie  equation  and  Humble  formula  under  effective  stress  variations. 
Relationships  were  derived  in  this  study  which  will  be  helpful  in  more 
accurately  predicting  reservoir  fluid  and,  therefore,  energy  recovery. 

Chapter  II  provides  a  review  of  pertinent  literature  concerning 
the  interrelationships  among  and  functional  properties  of  resistivity, 
porosity,  and  permeability.  Chapter  III  develops  the  theoretical  basis 
on  which  this  current  work  rests.  Chapter  IV  presents  a  description  of 
the  equipment  utilized  and  the  experimental  methods  employed.  The  re¬ 
sults  obtained  from  the  experiments  are  analyzed  in  Chapter  V.  Also, 
relationships  are  developed  relating  resistivity,  porosity  and  per- 


II.  REVIEW  OF  LITERATURE 

II. 1  Deformation  Studies  Relating  to  Porosity  and  Permeability 

Deformation  behavior  of  rocks  and  its  effects  on  porosity  and  per¬ 
meability  have  been  the  subject  of  many  studies.  A  comprehensive  review 
of  literature  on  the  subject  has  been  provided  by  Scorer  and  Miller2  and 
Friedman. 3  While  Scorer  and  Miller  discussed  compressibility,  porosity 
and  permeability  as  functions  of  pressure,  Friedman  concerned  his  review 
mainly  with  the  influence  of  porosity  and  permeability  on  the  mechanical 
behavior  of  rocks. 

Terzaghi4  provided  the  basic  concept  of  effective  stress  on  which 
most  work  in  soil  mechanics  and  rock  mechanics  is  based.  Biot5  generalized 
the  concept  to  the  three  dimensional  case  and  presented  the  classical 
theory  of  three  dimensional  consolidation  of  porous  media. 

Geertsma6  generalized  these  theories  to  Include  viscoelastic  de¬ 
formations  in  porous  rocks.  He  obtained  mathematical  relationships  for 
pore  and  bulk  volumetric  changes  in  terms  of  the  effective  pressure.  He 
found  that  three  elastic  constants  and  three  viscous  constants,  besides 
porosity,  were  needed  to  describe  the  volumetric  variations  in  the  rock. 
His  experimental  results  showed  that  consolidated  sandstones  deformed 
essentially  elastically. 

Nur  and  Byerlee7  derived  an  exact  expression  for  an  effective 
stress  law  for  elastic  deformation  of  aggregates  with  pore  pressure. 

Their  work  demonstrated  the  exactness  and  usefulness  of  the  effective 
stress  concept. 
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Hall8  demonstrated  experimentally  that  rock  compressibility  in 
undersaturated  reservoirs  was  a  strong  function  of  overburden  load.  He 


showed  that  porosity  reductions  with  increasing  overburden  pressure  could 
significantly  effect  fluid*- in-place  calculations. 

Fatt9  obtained  similar  results  in  his  studies  on  porosity  and 
permeability  changes  in  sandstones  as  functions  of  overburden  pressure. 

He  concluded,  however,  that  porosity  changes  were  negligible,  but  per¬ 
meability  reductions  were  large,  enough  to  warrant  consideration.  He 
used  the  term  "effective  permeability"  to  represent  the  permeability 
downhole  at  high  pressure  conditions. 

Gangi10  determined  theoretical  variations  if  porosity  and  per¬ 
meability  with  applied  pressure  using  phenomenological  models  based  on 
the  Hertz  Theory  of  deformation  of  spheres.  He  obtained  nonlinear  ex¬ 
pressions  for  porosity  and  permeability  as  functions  of  the  externally 
applied  pressure  which  were  shown  by  him  to  exhibit  most  of  the  character¬ 
istics  of  experimental  data  obtained  by  others. 

Yoelin1 1  found  that  variations  in  bulk  compressibility  with 
pressure  were  usually  more  significant  than  changes  in  porosity.  He  dis¬ 
covered,  however,  that  some  highly  compressible  sandstones  reflected 
disproportionally  greater  pore  volume  variations  with  increasing  over¬ 
burden  pressure. 

Gray  et  al.12  in  their  work  with  rock  samples  from  the  geo- 
pressured  provinces  of  south  Texas,  observed  small  (1-2%)  pore  volume 
reductions  and  significant  (10%)  permeability  declines  over  a  7,000  psi 
effective  stress  range.  They  found  that  dynamically  determined  values 
of  porosity  bear  very  little  relation  to  porosities  determined  statically. 


Even  though  results  from  the  two  methods  were  closer  at  high  effective 
stresses,  considerable  discrepancy  still  existed.  They  also  found  that 
compaction  coefficients  were  much  lower  than  previously  assumed  and  that 
in  many  cases  the  rocks  were  considerably  inelastic. 

II. 2  Resistivity  Studies 

Literature  concerning  resistivity  of  formations  has  long  been 
associated  with  electric  logging.  Resistivity  and  spontaneous  potential 
measurements  were  the  sole  purpose  of  early  logging  techniques.13  The 
applications  of  resistivity  data  since  then  have  grown  and  have  been 
studied  by  many  investigators. 

The  classical  works  of  Archie114  and  Winsauer  et  al.ls  pervade 
most  literature  on  formation  resistivity.  Archie  developed  the  concept 
of  formation  factor  which  is  defined  as  the  ratio  of  the  resistivity  of  a 
100%  saturated  rock  to  the  resisitivity  of  the  solute.  He  observed 
from  experimental  data  that  a  reasonably  accurate  relationship  existed 
between  the  formation  factor  and  the  rock  porosity  at  atmospheric  pres¬ 
sure.  This  relationship  is  expressed  as: 

F  -  <fTm  2.1 

where 

F  *  formation  factor 

p  ■  porosity  (decimal) 

m  -  slope  of  the  log-log  linear  equation 
He  found  that  for  loosely  or  partly  consolidated  sands  of  the  Gulf  Coast 


m  values  varied  from  1.3  to  2.0. 


Winsauer  et  al . 1 5  performed  a  comprehensive  study  on  the  resis¬ 
tivity  of  brine  saturated  sands  in  relation  to  pore  geometry.  They 
developed  an  expression  for  the  resistivity  as  a  function  of  the  tor¬ 
tuosity  of  pore  channels  in  a  porous  rock.  Through  experiments  they 
derived  the  folJ owing  relationship  between  the  formation  factor  and 
porosity: 

F  -  0.62  <jf  2* 1  5  2.2 

Guyod16  studied  the  resistance  of  permeable,  brine  saturated  for¬ 
mations  from  a  more  academic  point  of  view.  He  considered  the  conductivity 
of  formations  dependent  on  several  factors:  1)  the  amount  of  conductive 
solids  in  the  rock;  2)  the  type  and  concentration  of  the  solute;  3)  the 
temperature  and  pressure  of  the  solute;  4)  the  geometry  of  the  pore 
space;  5)  the  shale  content.  Based  upon  his  observations  that  the  m 
exponent  in  Archie's  equation  (2.1)  tended  to  increase  from  1.3  to  2.2 
when  traversing  from  noncemented  to  highly  cemented  sandstones,  he 
labelled  the  exponent  the  "cementation  factor,"  a  term  which  is  currently 
used  in  literature. 

Wyllie  and  Rose17  examined  the  theoretical  basis  of  electric  log 
interpretation  and  attempted  to  quantify  the  cementation  factor  and 
porosity  parameters.  They  found  the  practical  range  of  m  varied  from  1.3 
to  3.0  but  emphasized  that  m  was  really  only  of  qualitative  value.  They 
realized  that  porosity  could  be  computed  only  if  m  was  specified,  and  con¬ 
cluded  that  tortuosity  data  obtainable  from  core  analysis  was  necessary. 

Owen18  agreed  with  Wyllie  and  Rose  and  stated  that  cementation 
factors  were  difficult  to  determine  because  of  the  complexity  of  the  re- 
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lationship  between  resistivity  and  pore  geometry.  He  used  a  synthetic 
pore  system  model  to  determine  m  values  and  found  that  they  corresponded 
with  the  observations  made  by  Archie  and  Guyod.  He,  like  previous  in¬ 
vestigators,  conducted  his  experiments  at  atmospheric  pressure. 

Wyble19  dealt  with  resisitivty  as  a  function  of  pressure.  In  his 
experimental  work,  he  raised  the  confining  pressure  on  sandstone  samples 
and  measured  the  resistivity  and  porosity.  He  obtained  values  for  the 
cementation  factor  and  Winsauer  coefficient  but  did  not  treat  them  as 
variables. 

Fatt20  also  viewed  resistivity  and  porosity  as  functions  of  pres¬ 
sure.  He  showed,  however,  that  the  cementation  factor  was  not  a  constant 
and  that  it  changed  with  confining  pressure. 

Glanville21,  in  his  work  on  formation  factors,  noted  increases  in 
the  value  of  m  as  the  effective  stress  on  the  rock  samples  was  increased. 
He  observed  these  variations  at  an  elevated  pore  pressure  and  with  chang¬ 
ing  confining  pressure. 

Redmond22  also  found  cementation  factor  values  to  increase  with 
pressure.  For  overburden  loads  in  the  range  of  10,000  to  15,000  psi,  he 
noted  that  cementation  factors  for  some  sandstones  were  as  high  as  8.3. 

He  also  observed  porosity  changes  from  7%  to  20%  over  a  16,000  psi  over¬ 
burden  pressure  range,  a  variation  which  was  higher  than  what  was  ob¬ 
served  by  other  investigators.  Finally,  Redmond  found  that  resistivity 
changes  with  pressure  were  less  pronounced  than  corresponding  changes  in 
permeability. 

Friedman,3  in  his  review,  addressed  the  relationship  between  re¬ 
sistivity  and  permeability.  He  stated  that  the  association  between 
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conductivity,  i.e.  the  reciprocal  of  resistivity,  and  the  permeability 
(or  hydraulic  conductivity)  was  natural  since  both  were  analogously  re¬ 
lated  to  the  ability  of  a  rock  to  permit  flow.  Attempts  to  show  such  an 
analogous  relationship,  however,  have  not  been  successful. 

Wyllie  and  Spangler23  reported  on  permeability  and  resistivity 
tests  made  on  Pennsylvanian  sandstone.  They  found  that  for  this  clean 
sandstone  resistivity  decreased  with  increasing  permeability.  They  con¬ 
cluded,  however,  that  this  relationship  reflecting  the  permeability- 
resistivity  analog  might  not  be  so  apparent  for  rocks  containing  shale 
inclusions,  mica  flakes,  or  secondary  cementation  effects. 

Boggus24  found  a  similar  basic  relationship  between  permeability 
and  resistivity  for  Berea,  Pecos,  and  Banders  andstones.  However,  he, 
like  Wyllie  and  Spangler,  conducted  his  tests  at  atmospheric  pressure 
under  zero  effective  stress  conditions. 

Tixier25  developed  a  method  for  evaluating  permeability  from 
electric  log  resistivity  gradients.  Although  charts  stemming  from  his 
studies  are  in  use,  they  are  applicable  only  for  hydrocarbon  reservoirs 
and  do  not  attempt  to  support  a  universal  analog  between  resistivity  and 
permeability.  Field  correlations  by  Lebourg  et  al.26  are  similarly 
applicable. 

Tellinghuissen2 7  presented  a  laboratory  method  for  simultaneous 
determination  of  rock  resistivity  and  permeability  at  variable  confining 
pressures.  His  results  merely  agreed  with  the  findings  of  other  re¬ 
searchers  who  studied  these  two  parameters  separately. 
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III.  THEORETICAL  CONSIDERATIONS 


III.l  Effective  Stress  Concepts 

The  computations  of  the  deformational  characteristics  of  rocks  as 
reflected  in  resistivity,  porosity,  and  permeability  changes  in  a  pro¬ 
ducing  reservoir  were  based  on  the  theories  developed  by  Terzaghi1*  and 
Biot5.  A  brief  explanation  of  the  theory  used  is  presented  here. 

A  change  in  the  pore  pressure  within  the  reservoir  produces  the 
same  effects  in  the  reservoir  as  an  increased  loading  would.  The  concept 
of  effective  stress  says  that  there  is  a  one  to  one  relationship  between 
external  load  and  pore  fluid  pressure.  Effective  stress  in  general  is 
defined  as: 


PF  -  P  -  P„  3.1 

L  O 

where 

PE  “  effective  stress 

P  ”  externally  applied  stress  (compression  positive) 
PQ  -  pore  pressure 

Although  the  form  of  the  effective  stress  equation  has  found 
widespread  acceptance  and  application,  further  studies  have  suggested 
that  the  relationship  should  be  as  presented  by  Biot: 

p  -  p  -  n  p„  3.2 

E  O 

where  ri  ■  constant  <  1 
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The  absolute  value  of  q  has  not  been  agreed  upon, 
however,  proposed  the  following  form  for  r|: 

n-x-i 

n 

where 

CR  »  rock  matrix  compressibility 
Cg  *  rock  bulk  compressibility 

For  most  natural  aggregates  CR<  <  CR  and  n  -  1.  Therefore,  equation  3.2 
reduces  to  3.1.  It  is  clear  from  either  equation,  however,  that  a  de¬ 
crease  in  the  pore  pressure  serves  to  increase  the  effective  compressive 
stresses,  thereby  leading  to  compaction.  Compaction  of  the  rock  in  turn 
leads  to  porosity  and  permeability  changes. 

In  a  reservoir,  this  effective  stress  concept  is  useful  in  re¬ 
presenting  the  rock  stresses  under  insitu  conditions.  Provided  lateral 
dimensions  of  the  reservoir  are  large  compared  to  its  thickness,  re¬ 
servoirs  are  assumed  to  deform  predominantly  in  the  vertical  plane.  For¬ 
mation  compaction  is,  therefore,  characterized  by  vertical  deformation  of 
the  reservoir  as  a  result  of  the  effective  stress  change  in  the  formation. 

As  a  result  of  fluid  withdrawal  resulting  in  formation  compaction, 
the  resistivity,  porosity,  and  permeability  vary.  These  rock  parameters 

can  be  related  directly  to  the  effective  stress  as  shown  in  Sections  HI. 2 
and  III. 3. 

III. 2  Resistivity  and  Porosity  Relationships 

Resistivity  as  a  useful  rock  parameter  was  utilized  by  Archie14  in 
his  development  of  the  formation  factor  as 
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where 

F  -  formation  factor 

•  resistivity  of  the  100%  brine  saturated  rock 
(ohm  meters) 

Rg  “  resisitivity  of  the  saturating  brine  (ohm  meters) 


Archie  and  later  on,  several  other  investigators15’16*17  found 
experimentally  that  the  formation  factor  could  be  related  to  the  porosity 
by 


F  -  a  <jTm  3.5 

where 

< P  ■  porosity  (decimal) 
a  =*  constant 
m  »  cementation  factor 


Equation  3.5  is  merely  a  generalized  form  of  equations  2.1  and  2.2. 

Values  for  the  coefficient,  a,  and  the  exponent,  m,  vary  depending  on  the 
rock  type.  Although  in  Archie's  equation  a  ■  1,  Winsauer  et  al.15 
specified  a  -  0.62  for  consolidated  sandstones.  Wyllie  and  Rose17  showed 
that  although  values  of  m  could  range  from  1  to  infinity,  the  variations 
encountered  in  practice  lay  well  within  the  limits  of  1.3  and  3.0. 

Both  Archie  and  Winsauer  obtained  their  results  on  rock  samples 
under  atmospheric  pressure  conditions.  Many  investigators9 ' 1 °* 1 9’ 22’ 27, 
however,  have  shown  that  permeability  and  porosity  are  functions  of 
pressure.  The  resistivity  and,  therefore,  the  formation  factor  in  turn 


i 


* 


have  also  shown  to  be  pressure  dependent. 2 0  *  2 1 ' 2 2 ’ 2 7  Two  of  these  same 
studies20’21  have  Indicated  that  the  cementation  factor,  m,  is  sensitive 


to  pressure. 


In  view  of  these  investigations,  Archie's  equation  can  be  modified 


as  follows: 


F(Pe)  -  a  ♦(Pe)"b(PE) 


For  the  equation  to  be  theoretically  correct,  boundary  conditions  must 
hold.  Therefore,  since  the  formation  factor  should  equal  1  for  100% 
formation  porosity,  a  ■  1.  Equation  3.6,  then,  can  be  put  into  the 
following  form: 


F(Pe)  -  0(PE)“m<PE) 


III. 3  Resistivity  and  Permeability  Relationships 

The  flow  of  current  and  the  flow  of  fluid  through  rocks  have 
often  been  termed  analogous.  Both  resistivity  and  permeability  are 
essentially  dependent  upon  the  same  flow  paths  and  subject  to  the  same 
geometric  considerations.  Therefore,  in  rocks  composed  predominately 
of  nonconductive  solids,  such  as  clean  sandstones  which  exist  in  the 
producing  horizons  of  a  geopressured  reservoir,  the  current  flow  should 
follow  the  same  tortuous  route  as  the  hydrodynamic  brine. 

Morita28  showed  that  in  the  laboratory  radial  flow  through  a 
cylindrical  rock  core  could  be  represented  by 


k  ■  ('ap)(lho  (r,Hl000) 

x  e 


where 


. 
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k  ■  permeability  of  the  rock  sample  (millidarcies) 

q  *  fluid  flow  rate  (cc/sec) 

U  ■  fluid  viscosity  (cp) 

L  »  sample  axial  length  (cm) 

AP  -  differential  pressure  across  core  (atm) 

Rj  “  correction  factor  for  decreased  flow  rate 
(dimensionless);  (“0.94  for  the  work) 

K 

e 

— »  geometric  correction  factor  (dimentionless); 
e'  (-1.00  for  this  work) 


As  an  analog  to  Morita's  whole  core  type  permeameter,  Telling- 
huissen27  used  the  following  relationship  for  computing  resistivity: 


Rj  "  (AV)(L)(Rf)  (Kfit)  (100> 


where 

R^  “  rock  resistivity  (ohm-meters) 
i  ■  current  flow  (amperes) 

AV  »  voltage  drop  radically  across  core  (volts) 

Appendix  I  discusses  the  reliability  of  resistivity  values  obtained  using 
this  equation. 

The  interrelationship  between  permeability  and  resistivity  has 
also  been  the  focus  of  several  investigations. 1 7,2  3* 2 4’ 2  5>  26  These 
studies,  however,  were  not  accomplished  under  varying  effective  stress 
conditions  and  no  attempt  was  made  to  obtain  a  relationship  between  the 
permeability  or  resistivity  and  the  pressure. 


-  - 
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The  interrelationship  between  permeability  and  resistivity 
can  be  approached  using  the  fluid  flow  and  current  analogy  discussed 
above  with  pressure  as  an  lndapendent  variable.  Permeability  has  often 
been  labelled  the  hydraulic  conductivity ,  a  term  now  adopted  by  the  U.  S. 
Geological  Survey.  Likewise,  resistivity  can  be  called  reciprocal 
electrical  conductivity.  Since  the  formation  factor  is  a  direct  function 
of  the  resistivity,  F  must  also  be  inversely  related  to  the  electrical 
conductivity.  Also,  since  both  permeability  and  formation  factor  are 
pressure  dependent  as  discussed  in  Section  III. 2,  the  analogy  between  the 
two  parameters  would  imply  that  a  rock  experiencing  a  decrease  in  per¬ 
meability  with  effective  stress  would  also  experience  a  corresponding  in¬ 
crease  in  formation  factor.  Intuitively  then,  the  F  x  k  product  for  a 
particular  clean  sandstone  should  be  constant  under  changing  effective 
stress  conditions.  This  concept  can  be  expressed  as 

F  x  k  -  C  3.10 

where 

C  •  constant  (dependent  upon  the  rock) 

Substituting  equation  3.7  into  3.10  and  solving  for  k  results  in 


k 


3.11 


The  applicability  of  these  two  equations  in  relation  to  experimental 
results  is  discussed  in  Section  V. 3. 


V 


IV-  EQUIPMENT  AND  EXPERIMENTAL  PROCEDURE 

IV. 1  Simultaneous  Property  System  (SPS) 

The  basic  triaxial  stress  loading  vessel  and  pressure  control 
system,  called  SPS,  which  was  developed  at  the  Center  for  Earth  Sciences 
and  Engineering  and  is  described  in  detail  by  Evans29,  was  used  for  all 
experiments  (Figure  3).  This  system  is  designed  to  make  measurements 
on  samples  under  confining  pressures  up  to  30,000  psi  and  has  Independent 
controls  for  axial  loads  up  to  60,000  psi  and  pore  pressures  up  to  30,000 
psi. 

The  equipment  has  the  capability  for  simultaneous  determination 
of  static  and  dynamic  moduli,  permeability,  porosity,  and  resistivity 
of  rock  samples  at  various  stress  levels.  Static  moduli  including 
bulk  compressibility,  uniaxial  compaction  coefficient,  matrix  com¬ 
pressibility,  Young's  modulus,  and  Poisson's  ratio  can  be  calculated 
from  radial  and  axial  strains  measured  by  two  rosette  strain  gauges 
cemented  onto  the  rock  sample  (Figure  4) .  Stresses  are  measured  by  BLH 
pressure  transducers.  Dynamic  parameters  can  be  calculated  from  acoustic 
travel  times  for  p  and  s  waves  through  the  sample.  The  waves  are 
generated  and  received  by  peizo-elastic  crystals  contained  in  the  trans¬ 
mitting  and  receiving  transducer  cells  located  in  the  pressure  vessel 
heads  (Figure  3) . 

Permeabilility  and  porosity  measurements  were  made  using  the  MTS 
Intensifier  and  an  HIP  screw  pump.  The  intensifier  provides  constant 
upstream  pressure  to  cause  fluid  flow  through  the  rock  sample  for  per- 
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FIGURE  3  THE  SIMULTANEOUS  PROPERTY  SYSTEM  (SPS) 
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meability  determination  at  any  pore  pressure.  The  screw  pump  serves  to 
measure  pore  volume  changes  for  porosity  determination  under  elevated 
pore  pressures.  It  is  rated  at  30,000  psi,  has  a  capacity  of  11 cc,  and 
can  detect  pore  volume  changes  up  to  0.005  cc.  At  atmospheric  pore 
pressures,  volume  changes  are  determined  by  measurement  of  expelled  fluid 
into  an  inverted  pipette. 

Resistivities  were  calculated  from  measured  voltage  drops  across 
a  rock  sample  incurred  by  current  generated  with  a  General  Radio,  1310-B 
Oscillator.  Electrical  isolation  of  the  rock  sample  was  achieved  by 
epoxy  coatings  on  the  sample  ends,  by  delrin  insulation  fittings  around 
the  1/8  inch  high  pressure  tubing  pore  fluid  lines,  and  by  a  teflon 
packed  and  seated  high  pressure  valve  installed  in  the  external  by-pass 
line.  The  actual  resistivity  measurements  were  made  by  closing  the  by¬ 
pass  valve  and  observing  the  total  voltage  drop  in  the  resistivity  circuit 
and  the  voltage  drop  radially  across  the  cylindrical  rock  sample.  More 
detail  on  the  methods  employed  is  given  in  Appendix  I. 

IV. 2  Pore  Fluid  System 

The  pore  fluid  system  is  an  integral  part  of  the  SPS.  The  system 
can  be  divided  into  three  components:  1)  the  pumping  and  pressure  con¬ 
trol  equipment  needed  for  pore  pressure  attainment;  2)  the  core  jacket 
needed  for  containment  of  the  sample;  3)  the  pore  fluid  measuring  equip¬ 
ment  needed  for  permeability  and  porosity  determination. 

The  pressure  component  is  comprised  of  a  Sprague  pump,  HIP  screw 
pump,  MTS  intensifier,  BHL  pressure  and  Validyne  variable  reluctance 
transducers,  and  associated  Autoclave  valves  and  1/8  inch  stainless  steel 
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tubing.  The  detailed  schematic  of  the  system  is  available  in  Evans’  * 9 
original  work.  Calibration  data  for  the  apparatus  is  given  in  Appendix 


The  core  jacket  is  a  polyurethane  sleeve  equipped  with  two  pass¬ 
through  "shoes."  These  industrial  chromium  plated  stainless  steel  "shoes" 
are  molded  into  the  sleeve  and  serve  as  pore  fluid  passthroughs  from  the 
upstream  fluid  line  to  the  core  and  from  the  core  to  the  downstream  fluid 
line.  The  "shoes"  also  act  as  current  and  voltage  electrodes  for  re¬ 
sistivity  tests.  Figures  5  and  6  depict  the  core  jacket  and  passthrough 
arrangement  respectively.  Tellinghuissen2 7  provides  detail  on  the  basic 
jacket  and  "shoe"  design. 

The  pore  fluid  measurement  component  is  a  simple  arrangement  for 
determination  of  flow  rates  and  pore  volume  changes.  A  pipette,  manual 
timer,  and  the  MTS  intensifier,  which  maintained  a  constant  upstream 
pressure,  provide  necessary  data  for  calculation  of  permeabilities.  A 
system  of  pipettes  connected  to  the  downstream  lines  allow  the  measure¬ 
ment  of  pore  volume  changes  at  atmospheric  pore  pressure.  The  calibrated 
HIP  screw  pump  described  above  is  used  to  measure  pore  volume  changes  at 
elevated  pore  pressures. 


IV. 3  Sample  Preparation 

Sections  of  clean  sandstone  from  reservoir  depths  of  11,700  to 
15,700  feet  were  selected  from  4  inch  well  cores  obtained  from  Pleasant 
Bayou  wells  //I  and  if 2 .  Most  of  the  sections  tested,  however,  represented 
various  producing  zones.  These  sections  were  themselves  cored  perpen¬ 
dicular  to  the  apparent  bedding  planes  with  a  2  inch  I.D.  diamond  core 
barrel.  The  2  inch  cores  were  rough  cut  to  a  desired  length  with  a 
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Felker  diamond  cut-off  wheel  and  finally  cut  exactly  to  2.75  inches  on  a 
lathe  ensuring  that  the  end  faces  were  square  with  the  sides.  This 
1.375  length-to-diameter  ratio  was  selected  to  minimize  end  effects  in 
the  measurements . 2  9 

The  ends  of  the  samples  were  epoxied  with  Armstrong  A-12  ad¬ 
hesive  and  trimmed  on  a  lathe  to  a  0.004  inch  thickness  on  each  end. 

This  procedure  isolated  the  rock  sample  electrically  from  the  conductive 
heads  of  the  pressure  vessel  for  proper  resistivity  measurements.  The 
epoxy  did  not  penetrate  significantly  (which  was  checked  by  cutting 
coated  samples)  into  the  rock  to  affect  any  static  measurements.  How¬ 
ever,  its  effect  on  dynamic  p-wave  travel  times  was  corrected  for  in  the 
final  computations  of  p-wave  velocity  through  the  sample. 

In  order  to  obtain  bulk  volumetric  changes  of  the  samples  with 
pressure,  two  2-element  90°  'tee’  rosette  strain  gauges  were  cemented 
diametrically  opposite  each  other  on  the  rock  circumference.  The 
rosettes  were  Micro-Measurement,  1/8  inch,  350  ohm  gauges.  The  varnish 
insulated  lead  wires  were  covered  with  plastic  tape  to  prevent  any 
possible  shorting  due  to  brine  penetrating  the  insulation  (Figure  4) . 

After  the  gauging  process,  in  order  to  saturate  the  rocks,  the 
samples  were  subjected  to  a  vacuum  for  two  to  four  hours  depending  on 
their  Individual  permeabilities.  The  vacuum  process  was  followed  by 
flooding  the  samples  with  6%  NaCl  solution  which  ensured  that  the  samples 
were  saturated  fully  and  free  of  air  bubbles.  The  concentration,  as 
stated  by  Ucok30,  is  in  the  range  of  geothermal  brines.  It  also  gave 
pore  fluid  resistivities  at  room  temperature  which  were  low  enough  to 
allow  easy  detection  of  changes  in  the  rock  resistivity.  The  6 X  con- 


centration  also  was  found  to  be  helpful  in  preventing  clay  swelling 
within  samples. 
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When  ready  for  testing,  a  prepared  sample  was  measured  with  dial 
calipers  and  inserted  into  the  core  jacket.  The  sample  and  core  jacket 
were  then  clamped  between  the  axial  heads  of  the  pressure  vessel.  The 
strain  gauges  lead  wires  were  connected  to  panels  leading  to  the  Wheat¬ 
stone  bridge  circuitry  of  two  Bell  and  Howell  Mode  8-114  Signal 
Conditioners.  Upstream  and  downstream  pore  fluid  lines  were  fitted  to 
the  passthrough  shoes  of  the  core  jacket  and  leads  for  radial  resistivity 
measurements  were  attached  to  the  fittings  with  alligator  clips.  Finally, 
the  vessel  heads  and  sample  were  lowered  into  the  vessel  and  secured. 
Pressure  loading  and  rock  testing  followed. 

Initial  porosity  measurements  were  made  on  smaller,  1  inch  dia¬ 
meter  samples  cut  adjacent  to  the  respective  test  samples.  These  1  inch 
samples  were  oven  dried,  weighed,  vacuumed,  saturated,  and  lastly  re¬ 
weighed  to  obtain  porosity  values  at  atmospheric  conditions.  The  initial 
porosity  was  used  in  calculations  as  the  reference  from  which  variations 
in  this  parameter  were  measured. 


IV. 4  Test  Procedure 

Each  of  the  twelve  rock  samples  were  subjected  to  two  different 
pressure  loadings:  1)  hydrostatic  loading  in  which  pore  pressure,  con¬ 
fining  pressure,  and  axial  pressure  were  equal  and  increased  or  decreased 
simultaneously;  2)  trlaxial  loading  in  which  pore  pressure,  confining 
pressure,  and  axial  pressure  were  not  equal.  The  hydrostatic  loading  was 
used  basically  to  put  the  sample  initially  at  an  elevated  stress  level 
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from  which  different  effective  stress  conditions  (triaxial  loadings) 
were  generated.  Three  types  of  triaxial  loading  were  run  on  the  sample: 

1)  hydrostatic  compression  test  at  elevated  pore  pressure;  2)  uniaxial 
compaction  test;  3)  hydrostatic  compression  test  at  atmospheric  pore 
pressure  (not  preceeded  by  a  hydrostatic  loading) .  These  tests  are 
described  below. 

The  procedure  for  the  hydrostatic  compression  test  at  elevated 
pore  pressure  was  to  raise  the  pore,  confining,  and  axial  pressures  to  a 
value  corresponding  to  the  sample  depth  using  a  1  psi/ft  overburden 
gradient.  Then,  as  shown  in  Figure  7,  with  the  confining  and  axial 
pressures  kept  constant,  the  pore  fluid  pressure  was  reduced  in  steps 
to  one-half  the  specimen  depth  and  then  increased  back  to  the  initial 
value.  In  order  to  check  the  effective  stress  concept,  a  similar  variable 
pore  pressure  test  was  also  conducted  using  a  pressure  level  correspond¬ 
ing  to  psi  per  foot  times  the  sample  depth. 

For  the  uniaxial  compaction  test  (or  constant  radial  strain 
test),  the  sample  initially  was  hydrostatically  loaded  as  above.  Then 
as  depicted  in  Figure  8,  the  pore  pressure  was  reduced  in  steps  while 
keeping  the  axial  pressure  (overburden)  constant  and  while  adjusting  the 
confining  pressure  so  that  radial  (lateral)  strain  remained  constant.  As 
done  in  the  compression  test  above,  the  uniaxial  compaction  test  was  re¬ 
peated  using  a  pressure  level  corresponding  to  a  gradient  of  H  psi/ft. 

In  the  hydrostatic  compression  test  at  atmospheric  pore  pressure, 
no  initial  hydrostatic  loading  occurred.  As  shown  in  Figure  9,  the  con¬ 
fining  and  axial  pressures  were  raised  simultaneously  in  steps  to  the 
sample  depth  using  a  1  psi/ft  overburden  gradient.  The  pore  pressure 
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FIGURE  7  PRESSURE  LOADING  SEQUENCE  FOR  HYDROSTATIC  COMPRESSION 
TEST  AT  ELEVATED  PORE  PRESSURE 
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meanwhile  was  maintained  at  an  atmospheric  level.  When  a  pressure  cor¬ 
responding  to  the  sample  depth  was  achieved,  the  confining  and  axial 
pressures  were  reduced  in  steps  back  to  their  initial  values. 

The  three  trlaxial  loading  tests  described  above  were  used  to 
check  the  validity  of  the  effective  stress  concept  and  to  determine  if  sig¬ 
nificant  discrepancies  existed  among  porosity,  resistivity,  and  per¬ 
meability  values  when  they  were  obtained  under  different  loading  con¬ 
ditions.  These  parameters  were  measured  at  each  pressure  step  which  also 
permitted  their  dependence  on  effective  stress  to  be  observed. 

The  porosity  was  obtained  from  both  static  and  dynamic  measure¬ 
ments.  The  static  methods  utilized  volumetric  strains  and  pore  fluid 
volume  changes  for  calculation  of  porosity  values .  Another  static  method 
used  resistivity  measurements  which  allowed  porosity  calculations  with  the 
Humble  formula  (equation  2.2).  The  dynamic  method  utilized  observed  p- 
wave  travel  times  through  the  samples  which  could  be  used  in  the  Wyllie 
Time  Average  Equation  (equation  5.7)  to  compute  porosity  values. 

The  resisitivity  was  obtained  from  voltage  drops  measured 
radially  across  the  sample.  As  discussed  in  Appendix  I,  resistivities 
were  converted  to  formation  factors  using  pore  fluid  resistivities 
corrected  for  pressure.  The  formation  factors,  then,  were  employed  in 
data  analysis. 

The  permeability  was  calculated  using  observed  flow  rates  of  brine 
through  the  sample  under  a  differential  pressure  driving  force.  Equation 
3.8  for  a  whole  core  type  permeameter  was  utilized.28 


FIGURE  9  PRESSURE  LOADING  SEQUENCE  FOR  HYDROSTATIC  COMPRESSION  TEST 
AT  ATMOSPHERIC  PORE  PRESSURE 
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V.  ANALYSIS  OF  EXPERIMENTAL  RESULTS 


V. 1  Resistivity  Changes  with  Effective  Stress 

Resisitivity  changes  can  best  be  studied  in  terms  of  formation 
factor  variations.  The  formation  factor,  F,  in  contrast  to  the  resis¬ 
tivity,  is  independent  of  the  electrical  properties  of  the  saturating 
fluid.  Therefore,  changes  in  F  correspend  only  to  changes  in  the  in¬ 
trinsic  properties  of  the  rock  and,  as  explained  in  Chapter  III,  to  the 
pressure. 

From  resistivity  data,  formation  factors  were  computed  at  each 
pressure  step  in  the  triaxial  loading  tests  described  in  Section  IV. 4  and 
were  plotted  against  effective  stress.  The  validity  and  accuracy  of 
the  effective  stress  concept  was  checked  by  comparing  the  changes  in  F  in 
each  loading  test  for  the  same  effective  stress  change. 

Figures  10  through  13  show  comparisons  between  variations  of  F 
with  effective  stress  for  the  hydrostatic  compression  tests.  Some  dif¬ 
ference  in  curve  shapes  between  tests  run  at  atmospheric  pore  pressure 
and  tests  run  at  elevated  pore  pressures  can  be  observed.  This  difference 
was  probably  caused  by  two  phenomena:  1)  the  opening  of  microcracks  at 
elevated  pore  pressures;  2)  matrix  compressibility.  The  first  factor 
resulted  in  a  larger  and  more  uniform  change  in  F  when  the  microcracks 
closed  with  increasing  effective  stresses.  The  second  factor,  matrix 
compressibility,  was  accentuated  at  elevated  pore  pressures.  With  a 
larger  matrix  compressibility,  n  in  equation  3.3  decreased,  and,  there¬ 
fore,  the  actual  effective  stress  on  the  rock  sample  was  reduced.  Again 


Pa  =  Pc=  7.35  kpsi 


FIGURE  10  COMPARISONS  OF  VARIOUS  HYDROSTATIC  COMPRESSION  TESTS  (#14699) 
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IGURE  11  COMPARISONS  OF  VARIOUS  HYDROSTATIC  COMPRESSION  TESTS  (#15688) 


FIGURE  12  COMPARISONS  OF  VARIOUS  HYDROSTATIC  COMPRESSION  TESTS  (//14702 
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the  result  was  a  larger  more  uniform  change  in  F  for  tests  run  at  elevated 
pore  pressures. 

The  matrix  compressibility  and  the  microcracks  of  the  samples  do 
not  significantly  affect  variations  of  F  in  the  atmospheric  pore  pressure 
tests.  The  flatter  curves  at  the  higher  effective  stresses  support  this 
fact.  Therefore,  data  gathered  from  compression  tests  under  atmospheric 
pore  pressure  conditions  may  not  be  accurate  enough  for  predicting  rock 
properties  under  insitu  rock  conditions. 

Offsetting  among  the  curves  could  be  attributed  to  several 
reasons:  1)  Permanent  set  in  the  samples  caused  by  loading  and  unloading 
cycles  during  different  tests  on  the  sample;  2)  Temperature  change 
between  each  test  run.  (Appendix  I  discusses  that  an  average  room 
temperature  was  used  in  determining  the  resistivity  of  the  saturating 
fluid.) 

Another  comparison  can  be  made  between  the  hydrostatic  compression 
test  and  the  uniaxial  compaction  test.  Figures  14  through  18  show  typical 
plots  of  F  verses  effective  hydrostatic  stress,  for  compression  tests  at 
elevated  pore  pressures,  and  effective  axial  stress,  for  uniaxial  com¬ 
paction  tests.  The  good  agreement  among  the  curve  shapes  demonstrates 
that  the  formation  factor  is  similarly  influenced  in  these  two  type  tests 
by  changes  in  the  effective  stress.  The  pressure  level  did  not  appear  to 
be  a  strong  factor  influencing  the  data.  This  resulted  because  the  pore 
pressures  were  high  enough  initially  to  open  microcracks  and  to  compress 
the  rock  matrix.  Any  offsetting  of  the  curves  was  caused  by  the  factors 
described  above. 

Since  both  of  the  trlaxial  loading  tests  under  elevated  pore 
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pressure  conditions  appeared  to  give  comparable  results,  further 
analysis  with  other  samples  was  conducted.  Figures  19  through  24  show 
plots  of  F  verses  effective  stress  for  the  uniaxial  compaction  tests 
for  many  of  the  rocks  tested.  The  graphs  for  each  rock  show  that  the 
triaxial  loading  tests  at  various  elevated  pore  pressures  give  similar 
results.  The  near  closeness  of  the  curve  shapes  demonstrates  the  validity 
and  accuracy  of  the  effective  stress  concept. 

The  data  obtained  in  the  experiments  showed  that  formation 
factors  resulting  from  effective  stress  changes  are  nonlinear  and  could 
be  represented  by  a  polynomial  expression  in  pressure.  Curvilinear 
bivariate  regression  on  the  formation  factor  and  effective  stress  data 
for  each  test  was  performed.  A  simple  quadratic  expression  in  pressure 
of  the  form 

F  -  m  +  n2P£  +  n3Pg2  5.1 

where 

n^  -  coefficients  dependent  upon  the  rock  type 

was  tried  but  did  not  fit  the  data.  A  very  good  correlation  for  all 
rock  samples  was  obtained  by  regressing  In  F  and  Pg,  resulting  in  a 
second  degree  polynomial  of  the  form 

In  F  ■  m  +  n2Pg  +  Q3PE2  5.2 

Figure  25  shows  a  typical  curve  for  one  of  the  rock  samples  and  de¬ 
monstrates  the  goodness  of  fit  of  equation  5.2.  Data  curves  for  the 
remainder  of  the  samples  are  given  in  Figures  26  through  36. 
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FIGURE  16  COMPARISONS  OF  HYDROSTATIC  COMPRESSION  TEST  AND  UNIAXIAL 
COMPACTION  TEST  (//14703) 
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FIGURE  19  COMPARISONS  OF  VARIOUS  UNIAXIAL  COMPACTION  TESTS  (#14699) 
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FIGURE  21  COMPARISONS  OF  VARIOUS  UNIAXIAL  COMPACTION  TESTS  (#15668) 
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FIGURE  22  COMPARISONS  OF  VARIOUS  UNIAXIAL  COMPACTION  TESTS  (014711) 
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FIGURE  23  COMPARISONS  OF  VARIOUS  UNIAXIAL  COMPACTION  TESTS  (#14712) 
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FIGURE  25  COMPARISON  OF  FORMATION  FACTORS  OBTAINED  FROM  MEASURED 
RESISTIVITY  DATA  AND  FROM  EQUATION  5.2  (0 14765,  UNI¬ 
AXIAL  COMPACTION  TEST) 
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The  n^  coefficients  are  dependent  upon  the  rock  matrix,  amount 
of  conductive  solids  in  the  matrix,  shale  content,  tortuosity,  and 
porosity.  The  coefficients  varied  from  rock  to  rock,  but  for  the  sand¬ 
stone  samples  representing  the  Pleasant  Bayou  reservoir,  they  had  a  low 
variation  from  the  mean  as  discussed  in  Section  V.4. 


V. 2  Porosity  Changes  with  Effective  Stress 

During  the  compaction  process,  as  discussed  earlier,  both  the 
formation  factor  and  the  porosity  change  with  the  effective  stress. 

In  order  to  develop  a  proper  relation  between  porosity  and  F,  porosity 
changes  were  measured  as  the  effective  stress  was  varied.  Porosity  was 
obtained  by  both  static  and  dynamic  methods  as  described  below. 

A.  Static  Measurements 

1.  Porosity  by  Water  Volume  and  Strain  Measurement  Method 
In  this  method,  expelled  water  volume  measurements  pro¬ 
vided  pore  volume  changes,  and  volumetric  strain  measurements  repre¬ 
sented  bulk  volume  variations.  The  measurements  were  used  to  calculate 
porosity  via  the  following  expression: 


4>  o  V  o  -  Av_ 


Vo 


5.3 


eff 


where 


<)>  i  »  measured  porosity  at  each  stress  level 
<Po  “  initial  porosity  of  sample 
Vo  ■  initial  bulk  volume 

AV  ■  expelled  water  volume  from  rock  pores  during  de¬ 
formation 


A 


eff 


volumetric  strain  due  to  effective  stress 
(“  S  +  2eR> 


* 
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e  «  axial  strain 
z 

eR  =  radial  strain 

2.  Porosity  by  Water  Volume  Measurement  Method 

Assuming  that  bulk  volume  changes  are  direct  reflections 
of  pore  volume  changes,6  porosity  can  be  expressed  by: 


<p2 


4>oVo  -  AVw 
Vo  -  AV 

v 


5.4 


3.  Porosity  by  Volumetric  Strain  Measurement  Method 

As  a  counterpart  to  method  2,  assuming  that  pore  volume 
changes  are  attributable  mainly  to  bulk  volume  changes,6  the  porosity 
can  be  calculated  by  using  only  measured  strains: 


413 


5.5 


4.  Porosity  by  Resistivity  Measurement  Method 

In  this  method,  the  Humble  formula  (equation  2.2),  which 
is  used  for  Gulf  Coast  sandstones,  was  utilized.  It  is  given  by: 


r  -  °-62 
F  “  X2TTT 


or  upon  rearrangement 


P4  "  e 


,0.62. 
1n__F__ 
2.15 


2.2 


5.6 


B.  Dynamic  Measurements 

Porosity  by  Willie  Time  Average  Method:37 
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$5 


5,7 


where 

v  ”  p-wave  velocity  through  the  saturating  fluid  (for 
the  6%  NaCl  solution  used,  v^  -  5195  ft/sec) 

=  measured  p-wave  velocity  through  the  sample 

=  p-wave  velocity  through  the  matrix  material  (for  the 
sandstone  samples  selected,  v  »  19,600  ft/sec 
based  on  values  obtained  in  t¥ie  laboratory  for 
fused  quartz.29) 

Porosities  obtained  using  these  methods  were  plotted  against 
the  applied  effective  stress  for  each  of  the  tests.  Figures  37  through 
48  show  such  plots  for  the  various  rock  samples  subjected  to  similar 
effective  stress  conditions.  Each  of  the  graphs  shows  that  the  per¬ 
centage  change  in  porosity  per  1000  psi  effective  stress  change  is  de¬ 
pendent  upon  the  method  of  determination.  Table  1  provides  a  relative 
comparison  among  the  methods. 

<f>5  gave  the  largest  percentage  change.  Variations  in  this 
method  appear  in  general  to  be  excessive.  The  greatest  change  occurred 
at  low  effective  stress,  however,  and  in  the  4,000  to  7,000  psi  range  of 
effective  stress,  the  variation  with  pressure  decreased,  providing 
values  of  porosity  which  for  most  rocks  approximated  porosities  deter¬ 
mined  by  other  methods. 

The  porosity  obtained  by  the  water  volume  and  volumetric  strain 
method,  <p i,  was  considered  to  be  the  most  accurate  measure  of  the  true 
porosity.  The  calculations,  which  were  corrected  for  brine  compressibility 
variations  according  to  Gibson,31  incorporated  no  assumptions  beside 
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TABLE  1 

Relative  Comparisons  of  Porosity  Reductions 
With  Effective  Stress 


Method  of  Porosity 
Determination 

Symbol 

Average  %  Change  in  Porosity 
per  kpsi  Effective  Stress 

Change 

Water  Volume  and 
Volumetric  Strain 

<J>1 

0.7  -  1.1 

Water  Volume 

4>2 

0.6  -  1.0 

■J 

Volumetric  Strain 

<P3 

0.1  -  0.3 

Resistivity 

4*4 

0.5  -  0.9 

Time  Average 

<p5 

3.0  -  4.0 

Porosity  range  for  the  samples  tested  was  16  to  21% 


'I  i  1/  jfa  »*<■ 
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rock  isotropy  and  homogeneity. 

In  contrast,  both  the  porosity  by  water  volume,  $2,  and  the 
porosity  by  volumetric  strain,  <j>3»  methods  were  computed  based  on  the 
simplifying  approximation6  that  if 

CR  <  <  S 

where 

CR  =  rock  matrix  compressibility 

Cg  *  rock  bulk  compressibility 

then,  pore  volume  changes  approach  bulk  volume  changes.  Or  more 
directly, 

AVp  -  AVB 

where 

AVp  =*  change  in  pore  volume 

AVg  »  change  in  bulk  volume 

The  porosity  as  determined  by  the  water  volume  method,  <p2,  has 
been  used  widely  in  literature6 ’ 1 5 ’ 2 0 ’ 2 1 ’ 22  to  represent  porosity 
variations  with  pressure.  From  a  macroscopic  viewpoint,  the  small 
difference  between  <j>2  and  $1  shown  in  Table  1  could  justify  the  use  of 
the  more  easily  obtainable  $2.  However,  this  difference  could  be  sig¬ 
nificant  when  modelling  reservoirs  of  high  initial  porosity  and  ab¬ 
normally  high  formation  pressures. 

The  volumetric  strain  porosity,  <f>3,  reflected  percentage  changes 
in  porosity  which,  in  relation  to  results  from  the  other  methods,  were 
low.  The  actual  variations  in  the  rock  porosity  were  masked  by  using  re- 


latively  small  percentage  changes  In  the  bulk  volume  to  reflect  pore 
volume  changes. 

Porosities  calculated  by  the  resistivity  method  using  the  Humble 
formula,  $4,  reflected  changes  with  effective  stress  which  were  less  than 
those  given  by  The  Humble  formula,  then,  did  not  appear  to  be  com¬ 

pletely  accurate  in  providing  porosities  from  resistivity  data  at  the  high 
effective  stresses  found  in  the  geopressured  reservoir. 

Logically,  since  both  resistivity  and  porosity  vary  with  pressure, 
the  cementation  factor  in  the  Humble  formula  and  in  Archie's  equation 
should  also  change.  As  explained  in  Section  111.2,  several  investi¬ 
gators20'21’22  found  that  the  cementation  factor  was  not  constant  with 
pressure.  In  view  of  this  fact,  a  modified  Archie's  equation  may  be 
written  as 


F(Pe)  -  (KPE)“m(PE) 

3.7 

or 

lnF(P£)  -  -m(PE)ln<KPE) 

5.8 

Formation  factors  were  calculated  using  measured  values  of  rock 
resistivity  and  pressure-corrected  values  of  brine  resistivity  as  ex¬ 
plained  in  Appendix  I.  Porosities  were  determined  by  the  expelled 
water  volume  and  volumetric  strain  method  (method  1). 

Graphs  of  InF  versus  lnc|>i  were  nonlinear  which  reflected  the 
pressure  dependence  of  the  cementation  factor.  Therefore,  in  order  to 
evaluate  m(P^)  formation  factors  were  plotted  against  $1  for  the  dif¬ 
ferent  rock  samples  tested  under  various  loading  conditions.  Figure  49 
through  60  show  the  results. 
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The  data  was  best  fitted  by  a  second  degree  polynomial  in  P£  of 


the  form 


where 


-m(Pg)  »■  mi  +  nuPj.  +  mjPg 


m^  ■  constant  coefficients 


5.9 


Substitution  into  equation  3.7  results  in 


F(PE)  -  <KPE)niWE*n3PE2 


5.10 


Notice  that  for  P£  •  0,  Archie's  equation  is  reproduced. 

The  coefficients  in  the  quadratic  expression  of  the  cementation 
factor  are  a  function  of  the  rock  type.  Guyod16  explained  that  m  it¬ 
self  is  dependent  upon  the  degree  of  cementation  in  the  rock.  The 
coefficients,  m^,  are  similarly  dependent.  Their  values  fluctuated 
slightly  from  sample  to  sample  tested,  but  the  statistical  variation  from 
the  mean  was  low  as  shown  in  Table  2  and  discussed  in  Section  V.4 

An  expression  for  the  porosity  as  a  function  of  the  effective 
stress  was  formulated  using  the  mean  values  of  the  coefficients  in  the 
quadratic  expression  of  m.  From  equation  5.10 


lnF(PE) 

ln<f>(PE)  -  mi  +  m2PE  + 


Substituting  equation  5.2  and  solving  for  porosity  results  in 


ni  +  n2P_  +  naP, 


5.11 


»i  +  m2Pg  +  nuPg 


5.12 


14,702  II 
Pe  •  14  7  kpit 


72 


EFFECTIVE  AXIAL  STRESS  (kpii)  EFFECTIVE  AXIAL  STRESS  (kp«i) 


9911 
7.35  kpti 


EFFECTIVE  AXIAL  STRESS  Up«l) 


V 


EFFECTIVE  HYOROSTATIC  STRESS  (kpii) 


FIGURE  41A 


EFFECTIVE  AXIAL  STRESS  (kpti ) 


FIGURE  4 IB 


14,712ft 


14,696 A  ft 


EFFECTIVE  AXIAL  STRESS  ( kpsi ) 


FIGURE  46 


DEPTH 
#14,711  It 


EFFECTIVE  AXIAL  STRESS  (kpti)  EFFECTIVE  AXIAL  STRESS  (Mil) 


EFFECTIVE  HYDROSTATIC  STRESS  (kpti) 


FIGURE  48B 


FORMATION  FACTOR,  F  2  FORMATION  FACTOR 


VTVTTYl 


IGURE  51A  FIGURE  51B 


Z3 


1953“ 


POROSITY,  <f> 


FIGURE  SIC 


FIGURE  5 ID 


FORMATION  FACTOR 


FIGURE  53B 


FORMATION  FACTOR 


(JdVOOOO  ♦ 


FIGURE  59C  FIGURE  59D 


FORMATION  FACTOR.  F  2  FORMATION  FACTOR 


FIGURE  60B 


<► 


•V 


s 


94 


A  change  in  porosity,  d4>,  as  a  result  of  an  effective  stress 


change,  dPg,  therefore  is  given  by 

Cr 


d<p 


( n :  +  2n3P£) 


(mi  +  m2PE  +  msPj/)  “ 


ni  +  n2PE  +  mPE 
mi  +  m2P£  +  m3PE2 


(m2  +  2m3PE)(ni  +  n2Pg  +H3PE 


dP. 


(mi  +  m2P£  +  maPj,2)2 


5.13 


Results  produced  by  equation  5.13  were  comparable  to  experimentally 
obtained  data.  Figure  63  shows  calculated  and  measured  values  of  porosity 
for  a  typical  sample.  An  analysis  of  the  results  is  provided  in  Section 
V.4. 


V.3  Permeability  Changes  with  Effective  Stress 

The  permeability  changes  with  effective  stress  observed  in  this 
work  were  qualitatively  similar  to  those  found  by  other  inves¬ 
tigators9’  1  0 ’ 12 ’ 22 ’ 2 7.  Permeability,  k,  was  found  to  decrease  signi¬ 
ficantly  with  increases  in  effective  stress.  Figure  61  shows  such  a 
dependence  for  one  of  the  rocks  tested.  All  other  samples  gave  similar 
results. 

A  general  polynomial  expression  of  k  in  P£  was  not  feasible 
because  of  the  non-uniform  variations  of  the  permeability  with  effective 
stress.  A  more  practical  and  obtainable  relationship  was  found  when 
permeability  was  related  to  resistivity. 

Wyllie  and  Spangler23  and  Boggus24  have  shown  that  resistivity 
increased  with  decreasing  permeability  for  clean  sandstones.  Since  re¬ 
sistivity  is  the  reciprocal  of  electrical  conductivity  and  permeability 


1 


4* 


AMiWihirfe.’  t 


has  been  labelled  the  hydraulic  conductivity,  this  relationship  seems 
to  be  natural. 

Wyllie  and  Spangler,  however,  concluded  that  this  relationship 
might  not  be  so  apparent  for  rocks  containing  shale  inclusions,  mica 
flakes,  or  secondary  cementation  effects.  Rock  samples  selected  for  this 
work  from  the  producing  horizon  of  Pleasant  Bayou  wells  41  and  42,  were 
relatively  clean  sandstones  and,  therefore,  should  not  have  been  in¬ 
fluenced  greatly  by  these  factors. 

Data  was  analyzed  using  the  hydraulic  -  electrical  conductivity 

analog  developed  in  Chapter  III  and  reflected  in  equation  3.10: 

F  x  k  -  constant  3.10 

The  product  of  permeability  and  formation  factor  at  each  stress  level 
for  each  rock  sample  was  normalized  in  analysis  to  permit  comparisons 
between  various  samples.  Since  experimental  data  were  not  always 
available  at  the  zero  effective  stress  condition,  a  5,000  psl  effective 

stress  level  was  used  as  the  reference  for  normalization. 

Figure  62  shows  the  normalized  Fxk  product  plotted  against 
effective  stress  for  the  twelve  rock  samples  tested.  Results  were  re¬ 
gressed  linearly  and  are  represented  by  the  equation 

y  -  1.000  -  0 . 002x  5.14 

where 

v .  <y  x  k) 

y  <F  x  k)  .  , 

reference  level 

Equation  5.14  shows  that  the  product  Fxk  is  essentially  independent  of 
the  effective  stress,  thereby  indicating  that 


V 
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Fxk 


(F  x  k) 


reference  level 


5.15 


Since  (yxk)reference  level  ^or  eac^'  roc^  *8  a  constant,  equation  5.15 
reduces  to  the  form  of  equation  3.10.  That  is 


Fxk-  constant  3.10 

Since  InF  can  be  represented  by  a  quadratic  expression  in  P£ 
(equation  5.2),  the  permeability,  upon  rearrangement  of  equation  3.10, 
can  be  given  by 

Ink  -  C  -  (m  +  n2PE  +  nsPg2)  5.16 

*  C  “  n2Pfi  -  n3PE2 
where  C  &  C*  are  constants 


The  permeability,  therefore,  becomes  a  function  of  the  rock  type  and 
the  effective  stress.  The  constant,  C',  can  be  evaluated  using 
permeability  data  for  a  core  at  atmospheric  pressure  where  P_  -  0,  re- 
suiting  in 


C' 


[Ink] 


PE-0 


5.17 


Equation  3.10  can  also  be  used  to  obtain  a  relation  between 
permeability  and  porosity.  Substituting  equation  5.10  into  3.10  re¬ 
sults  in 

k  x  ^mi+1B2l,E+mj^E  -  constant  -  C  5.18 


where  the  constant,  C,  depends  on  the  rock  characteristics.  Therefore, 
permeability  is  a  function  of  the  porosity,  effective  stress,  and  rock  type. 


01234567765432  I 
EFFECTIVE  AXIAL  STRESS  (kpsi) 


FIGURE  61  GENERAL  PERMEABILITY  TREND  WITH  EFFECTIVE  STRESS 
(#14711) 
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Comparisons  between  permeabilities  obtained  from  actual  experimental 
results  and  those  predicted  by  equation  5.18  are  given  in  Table  6  and 
Figure  64. 

V.4  Discussion  of  Results 

Values  for  coefficients  n^  and  as  expressed  in  equations  5.2 
and  5.9  respectively  were  calculated  for  each  test  (the  loading  pro¬ 
cedures  for  which  are  described  in  Chapter  III)  on  all  rock  samples. 

Table  2  lists  the  results.  One  can  observe  that  when  P£  «  0,  the  ce¬ 
mentation  factor  expression  reduces  to  a  constant,  mi,  with  a  value  in 
the  range  1.7  -  2.1.  These  values  are  consistent  with  those  given  by 
Archie  and  other  investigators.16’17 

A  statistical  analysis  of  the  results  provided  a  basis  for 
comparison.  Since  the  porosity  range  among  the  samples  was  small 
(162  -  212),  the  arithmetic  mean,  standard  deviation,  and  variation 
were  computed  for  each  coefficient.  Because  different  stress  loadings 
were  imposed  in  the  various  tests,  results  were  grouped  accordingly,  and 
the  variations  from  the  mean  compared.  Table  3  depicts  these  variations 
for  the  data  in  total  and  for  the  individual  tests.  In  general,  the 
lowest  variations  for  the  six  coefficients  occurred  for  the  uniaxial 
compaction  test. 

Since  uniaxial  compaction  loading  can  be  considered  to  best  re¬ 
present  the  stress  conditions  in  a  depleting  reservoir,  the  mean  and 
actual  m^  and  n^  values  for  this  test  were  used  to  back-cslculate  re¬ 
servoir  parameters,  <J>  and  k.  Changes  in  porosity  were  calculated  via 
equation  5.13  and  compared  to  actual  experimental  data.  All  comparisons 
were  made  using  <f>i  data  since  porosity  determined  by  the  water  volume  and 


i 
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volumetric  strain  measurement  method  was  most  reflective  of  the  true 
rock  porosity  as  discussed  earlier.  Tables  4  and  5  and  Figure  63  show 
‘  excellent  agreement  between  calculated  and  measured  porosity  changes 
for  several  rock  samples  when  the  actual  m^  and  n^  values  were  used. 
Although  some  offsetting  occurred  when  mean  coefficient  values  were 
utilized,  porosity  changes  with  effective  stress  variations  were  ap¬ 
proximately  the  same. 

Permeabilities,  back-calculated  using  the  mean  and  actual  m^ 
and  values  in  equation  5.18,  provided  a  method  of  verifying  the 
permeability-porosity  relationship.  Comparisons  were  made  with  ex¬ 
perimentally  measured  permeabilities.  Results  given  in  Table  6  show 
good  agreement  between  measured  k's  and  k's  calculated  from  actual 
coefficient  values.  Some  discrepancy  occurs,  however,  when  mean  co¬ 
efficient  values  are  used.  Figure  64  clearly  depicts  these  results 
for  a  typical  sample  during  a  uniaxial  compaction  test. 

In  general,  a  modified  Archie's  equation  and  a  permeability- 
porosity  relationship  as  given  by  equations  5.13  and  5.18  have  the 
capability  to  estimate  both  porosity  and  permeability  variations  with 
effective  stress. 


d$ 


n2+2n3P£ 

7 

(mi+m2PE+m3PEz) 

(mz+2m3Pg) (ni+nzPg+naPg2) 
(mi+nuPj.+msPg*) 


ni+n2PE+n3PE 

emi-Hn2PE+m3PE 


dPE 


5.13 


k  x  4.mi'h“2PE+BjPE2 


-  C 


5.18 
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TABLE  2 

Coefficients  of  the  Cementation  Factor  and  Formation  Factor  Expressions 


-<V  - 

mi  +  a2PE 

*  «iPEJ 

U,  P(PE) 

■  ni  nzPg 

+  ojPe!*** 

Sample* 

Test  Type** 

•  1 

*3 

na 

n2 

ni 

15665 

P  -P -7.85 

-1.9014 

.0011 

.0002 

1.3069 

.0104 

-.0005 

’A"''85 

-1.9508 

.0030 

.0001 

1.3394 

.0056 

-.0001 

11755 

P  -P.-12.0 

-1.9601 

-.0090 

.0014 

1.3004 

.0239 

-.0019 

*^u?rm  6*° 

-1.8616 

-.0174 

.0020 

1.3031 

.0225 

-.0017 

P*-1$.0 

-1.9863 

-.0105 

.0008 

1.3279 

.0223 

-.0015 

P^-6.0 

-1.9732 

-.0089 

.0010 

1.3808 

.0156 

-.0010 

14702 

P  -Pr-14.7 

-2.0413 

-.0036 

.0008 

1.3735 

.0096 

-.0004 

P*-P>7.35 

-1.9939 

-.0047 

.0010 

1.3776 

.0133 

-.0010 

P*-7?35 

-2.0133 

-.0014 

.0001 

1.3919 

.0100 

-.0003 

P^-14.7 

-2.0645 

-.0040 

.0007 

1.3869 

.0107 

-.0005 

14699 

P  -14. 7 

-1.8598 

-.0056 

.0006 

1.3596 

.0129 

-.0005 

P?-  7-35 

-1.7614 

-.0126 

.0010 

1.3319 

.0173 

-.0007 

P>P  -14.7 

-1.8394 

-.0076 

.0010 

1.3450 

.0146 

-.0008 

7-35 

-1.8108 

-.0088 

.0014 

1.3647 

.0162 

-.0011 

146968 

P  -P.-7.35 

-1.8199 

-.0094 

.0013 

1.3506 

.0178 

-.0009 

P*-lS.7 

-1.9291 

-.0040 

.0008 

1.4217 

.0163 

-.0011 

15665 

P  -  7.85 

-2.0836 

-.0010 

.0003 

1.4266 

.0080 

-.0002 

P*-15.7 

-2.1300 

-.0080 

.0014 

1.4095 

.0103 

-.0005 

P*-P  -  7.85 

-2.0708 

-.0039 

.0009 

1.4181 

.0129 

-.0009 

F>C-15.7 

-2.1172 

-.0061 

.0011 

1.3941 

.0199 

-.0019 

14712 

P  -P  -7.35 

-1.9064 

-.0030 

.0010 

1.3699 

.0133 

-.0011 

?*-  5. 35 

-1.9106 

-.0013 

.0005 

1.3745 

.0077 

-.0005 

P^-14.7 

-1.9633 

-.0056 

.0010 

1.3665 

.0099 

-.0005 

14696A 

P  -  7.5 

-1.8167 

-.0099 

.0010 

1.3549 

.0195 

-.0010 

PV-15.0 

-2.0206 

-.0085 

.0012 

1.4525 

.0156 

-.0011 

P7-Pr-15.0 

-1.9034 

-.0085 

.0008 

.1.3748 

.0138 

-.0006 

r  -p^-  7.5 

*1  L 

- 

- 

1.3680 

.0177 

-.0009 

14751 

P  -P.-7.0 

-1.9262 

-.0022 

-.0005 

1.3601 

.0079 

-.0002 

-1.9607 

-.0053 

.0011 

1.3624 

.0094 

-.0005 

14765 

Pa-14.765 

-1.7346 

-.0229 

.0023 

1.2390 

.0235 

-.0017 

14711 

P.-P.-  7.35 

-1.8590 

.0033 

.0055 

1.2566 

.0007 

-.0006 

P*"PSl4.7 

-1.8460 

.0039 

4.40x10  * 

1.2366 

.0117 

-.0008 

P*-  5.35 

-1.8773 

-.0017 

.0006 

1.2321 

.0083 

-.0005 

P^-14.7 

-1.9665 

.0015 

.0001 

1.2502 

.0103 

-.0006 

14703 

P  -Pr«12.0 

-1.8920 

.0054 

-.0006 

1.3917 

.0104 

-.0003 

^•7 

-1.9317 

-.0046 

.0007 

1.3967 

.0124 

-.0007 

*  Samp Is  number  reflects  depth  (ft*)  from  which  sample  obtained. 

«  Teat  type  la  described  by  the  Initial  hydrostatic  or  axial  stress  level,  reflecting 

respectively  a  hydrostatic  compression  or  uniaxial  compaction  test  as  discussed 
In  Chapter  III. 

***  n  values  are  based  upon  a  logia  derivation.  To  obtain  values  for  use  in  Ln# 
calculations,  multiply  listed  values  by  2.3. 
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TABLE  3 

Statistical  Parameters  for  the  Coefficients  of  Che  Cemencation 
Factor  and  Formation  Factor  Expressions 


Group  Analysis  Statistic* 

mi 

■i 

m3 

ni 

nj 

n3 

All  Tests 

y 

-1.934 

-0.005 

0.001 

1.353 

0.014 

-0.001 

a 

0.096 

0.006 

0.001 

0.055 

0.  005 

0.001 

v  (Z) 

5 

104 

111 

4 

34 

63 

•  * 

Uniaxial  Compaction 

u 

-1.959 

-0.007 

0.001 

1.361 

0.014 

-0.001 

Tests  8  1.0  psi/ft 

a 

0.103 

0.006 

0.001 

0.067 

0.005 

0.0004 

T 

Ini  dal  Overburden 
Pressure 

V  (Z) 

5 

as 

60 

5 

33 

44 

Uniaxial  Compaction 

y 

-1.928 

-0.005 

0.001 

1.355 

0.011 

-0.001 

Tests  @0.5  psi/fc 

a 

0.093 

0.004 

0.0004 

0.051 

0.005 

0.0003 

Initial  Overburden 
Pressure 

v(Z) 

5 

83 

61 

4 

41 

57 

Hydrostatic  Com¬ 

y 

-1.943 

-0.004 

0.001 

1.345 

0.015 

-0.001 

pression  Tests 

a 

0.096 

0.006 

0.001 

0.003 

0.005 

0.001 

(Elavated  Pore 
Pressure)  @1.0 
psi/ft  Initial 
Stress 

V  (Z) 

5 

151 

98 

4 

32 

65 

Hydrostatic  Com¬ 

u 

-1.906 

-0.005 

0.001 

1.348 

0.014 

-0.0009 

pression  Tests 

a 

0.079 

0.006 

0.002 

0.044 

0.004 

0.0004 

(Elevated  Pore 

v  (Z) 

4 

117 

112 

3 

31 

43 

J- 

Pressure)  i?  0.5 
pal/ ft  Initial 

Stress 


Statistical  Parameters  Are: 


L  *  Adthmetlc  Mean 
0  ■  Standard  Deviation 
V  •  Coefficient  of  Variation 


z  x  100 
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Mi 
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TABLE  4 

Comparison  of  Measured  end  Calculated  Porosities 
For  Several  Rock  Samples 


Sample 

Test 

P 

0 

m 

♦c 

A0 

m 

A0 

c 

A0 

c 

«1  4  nt 

“l  4  Qi 

■i 4  Di 

4  n 

Actual 

Mean 

Actual 

Mean 

14702 

PA 

14.7 

14.7  *  7.0 

.1995 

.2002 

.1906 

-.0136 

-.0133 

-.0167 

14711 

PA 

14.7 

14.7  -  7.35 

.2180 

.2171 

.1902 

-.0123 

-.0224 

-.0159 

14703 

PA 

14.7 

14.7  -  7.35 

.1766 

.1769 

.1902 

-.0115 

-.0150 

-.0159 

14699 

PA 

14.7 

14.7  -  7.0 

.1720 

.1716 

.1906 

-.0145 

-.0142 

-.0167 

14696B 

PA 

14.7 

14.7  -  7.0 

.1697 

.1687 

.1906 

-.0107 

-.0204 

-.0167 

Avarage: 

.1872' 

.1869 

.1904 

-.0125 

-.0171 

-.0164 

Notes:  1.)  $  &  AO  ■*  Measured  Data  Values 

m  m 

2. )  0C  4  A$c  ■*  Calculated  Values  Using  Eqns.  5.12  &  5.13  Respectively 

3. )  Porosities  Are  Given  @  P-[ 

c.  max. 

4. )  Porosity  Change*  Are  Civen  For  APF 


105 


TABLE  5 

Comparison  of  Measured  and  Calculated  Porosities 
For  a  Typical  Rock  Sample  ( //14696B) 


Sample 

Test 

Po 

(kpsi) 

<p 

rm 

*c 

mi  &  ni 

*c 

mi  &  ni 

Actual 

Mean 

14696B 

PA  -  14.7 

A 

14.7 

.1817 

.1832 

.2019 

12 

.1762 

.1763 

.1967 

10 

.1727 

.1723 

.1935 

7 

.1675 

.1687 

.1906 

10 

.1724 

.1723 

.1935 

12 

.1757 

.1763 

.1967 

14.7 

.1808 

.1832 

.2019 

Notes:  1.)  <p  ■*  Measured  Data  Values 
m 

2.)  <j>c  -+  Calculated  Values 
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FIGURE  63  COMPARISONS  OF  POROSITIES  OBTAINED  FROM  MEASURED  DATA  &  FROM  EQUATION  5.12  (#14696B 
UNIAXIAL  COMPACTION  TEST) 


'.£<1,  "J3&.  *'*  *  jr* 


v. 


TABLE  6 

Comparison  of  Measured  and  Calculated  Permeabilities 
For  Several  Rode  Samples 


Sample 

Test 

Po 

k 

vn 

kc 

kc 

Value  o 

(kpsi) 

(ad) 

«i  6  nt 

■i  &  “i 

Constat 

C 

Actual 

Mean 

14703 

Pa-14.7 

14.7  *  7.35 

20 

20 

21 

550 

14711 

Pa-14.7 

14.7  ♦  8.0 

156 

134 

103 

2635 

14696B 

PA-  14.7 

14.7  -  7.0 

87 

87 

102 

2635 

14699 

P.  -  14.7 

A 

14.7  -  7.0 

85 

96 

99 

2550 

14699 

PA  -  7.35 

6 

116 

113 

107 

2550 

A 

4 

107 

106 

102 

2550 

2 

101 

100 

100 

2550 

0 

95 

97 

99 

2550 

2 

100 

100 

100 

2550 

4 

104 

106 

102 

2550 

6 

108 

113 

107 

2550 

Notes:  1.)  k  •*  Measured  Data  Values 
01 

2. )  It  Calculated  Values  using  Eqn.  5.18  and  Calculated 

c  Porosities  from  Table  5 

3. )  k's  are  given  @  PF| 
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FIGURE  64  COMPARISONS  OF  PERMEABILITIES  OBTAINED  FROM  MEASURED  DATA  AND 
FROM  EQUATION  5.18  (#14699,  Uniaxial  Compaction  Test) 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


VI. 1  Concluaions 

1.  The  effective  stress  concept  is  a  valid  and  very  good  ap¬ 
proximation  of  the  stress  conditions  in  a  porous  rock.  However,  effective 
stresses  under  atmospheric  pore  pressure  conditions  do  not  adequately 
reflect  insitu  rock  conditions  and  may  not  be  accurate  in  predicting 
actual  reservoir  parameter  changes  with  pressure  if  the  rock  matrix 
compressibility  is  not  taken  into  account. 

2.  Porosity  changes  determined  by  using  only  measured  variations 
in  the  pore  fluid  volume  adequately  approximate  actual  porosity  changes 
in  medium  and  low  porosity  rocks  under  normal  stress  conditions.  How¬ 
ever,  for  highly  porous  rocks  or  for  a  more  accurate  Indication  of  the 
true  porosity  in  medium  porous  (~  20%)  rocks,  porosity  calculated  with 
both  volumetric  strain  and  pore  fluid  volume  measurements  should  be  used. 

3.  Resistivity,  porosity,  and  permeability  are  functions  of 
the  effective  stress  and  vary  as  the  stress  conditions  change.  Re¬ 
sistivity  and  permeability  are  more  sensitive  to  this  change  than  porosity 
and  showed  a  10-15%  variation  over  a  7,000  psi  pressure  range.  Porosity 
changes  reflected  a  5-7%  variation  over  the  same  range  for  the  sandstones 
tested. 

4.  The  cementation  factor  is  a  function  of  pressure.  In  log¬ 
ging  interpretation,  therefore,  a  modified  form  of  the  Archie  equation, 
as  developed  in  this  research  for  the  clean  sandstone  rocks  tested. 


must  be  used. 
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5.  Electrical  and  hydraulic  conductivities  are  analogous  for 
clean,  nonconductive  sandstones.  The  Fxk  product  is  independent  of 
pressure  for  a  given  rock.  This  relation  was  found  to  be  useful  in 
predicting  permeability  changes  with  pressure.  In  general,  such  pre¬ 
dictions  can  be  made  only  if  the  constant  of  proportionality  (equation 
3.10)  between  F  and  1/k  can  be  evaluated. 

6.  For  the  geopressured  zones  of  Pleasant  Bayou  wells  ft  1  and  // 2, 
the  permeability-porosity  relationship  is  log  linear  (equation  5.18)  and 
may  be  useful  for  modelling  purposes. 

VI. 2  Recommendations 

Based  upon  the  experience  gained  in  the  present  study,  the  follow¬ 
ing  suggestions  for  possible  improvements  are  made: 

1.  The  pore  fluid  temperature  should  be  iso thermally  controlled 
in  order  to  prevent  discrepancies  in  the  formation  factor  calculations 
due  to  pore  fluid  resistivity  changes. 

2.  Rocks  with  a  wider  range  of  porosity  and  permeability  should 
be  tested  in  the  SPS  in  order  to  verify  the  qualitative  results  of  this 
work. 

3.  An  Improved  method  of  permeability  measurement  should  be 
developed  to  more  easily  evaluate  changes  in  the  permeability  with  ef¬ 
fective  stress. 

4.  High  temperatures  reflecting  lnsitu  thermal  conditions 
should  be  incorporated  into  the  SPS.  Additional  tests  on  rock  samples 
should  be  conducted  in  equipment  presently  under  development. 
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A.  Resistivity  Measurement  Method 

Resistivity  measurements  in  this  study  were  made  from  observed 
potential  drops  across  the  flow  heads  as  described  earlier.  These 
potential  drops  are  a  function  of  not  only  the  rock,  saturating  fluid, 
and  current  flow  but  also  any  contact  resistance  that  exists  between 
the  electrodes  and  the  rock  and  any  polarization  of  charges. 

In  order  to  obtain  an  accurate  measurement  of  the  rock  re¬ 
sistivity,  both  contact  resistance  and  charge  polarization  must  be  elim¬ 
inated.  Dunlap  et  al.32  used  a  four-electrode  type  circuit  to  remove 
the  phenomena. 

In  this  method  current  was  passed  axially  through  the  sample 
using  end  plates  as  current  electrodes.  Potential  drops  were  measured 
using  additional  probes  or  electrodes  dispersed  along  the  sample  length. 
Dunlap  measured  and  allowed  for  any  appreciable  contact  resistance  by 
Inserting  a  high  impedence  resistor  of  the  same  size  as  the  input  im- 
pedence  of  the  voltmeter  in  series  with  the  meter.  (Figure  Al-1  simply 
depicts  the  circuitry.)  If  the  contact  resistance  were  negligible,  the 
new  voltage  after  insertion  of  the  resistor  would  be  one  half  the  pre¬ 
viously  measured  value.  If  the  contact  resistance  were  significant,  the 
new  voltage  would  be  greater  than  one  half  the  previously  measured  value. 

Later  work  as  described  by  Rust33  used  a  two-electrode  method  in 
which  the  same  electrodes  were  utilized  for  both  current  flow  and 
potential  drop  measurements.  Contact  resistance  between  the  electrodes 


111 


RAM  SI 
YSTQi 


FIGURE  Al-1  RESISTIVITY  CIRCUITRY  DUG! 

VOLTAGE  DROPS  WITHIN  THE  S' 
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and  the  brine  saturated  rock  sample  was  negated  by  several  methods  with 
varying  degrees  of  success.  Rust  used  conducting  silver  paint  on  the 
electrodes  to  reduce  the  contact  resistance.  His  results  showed  that 
provided  proper  techniques  were  employed  either  the  two-electrode  or 
four-electrode  method  yielded  reliable  resistivity  values. 

Dunlap  noted  that  polarization  within  the  sample  could  affect 
the  accuracy  of  the  measured  voltage  drops.  Although  alternating  current 
and  the  'iour-electrode  method  eliminate  any  significant  polarization 
problems,  AC  flow  introduces  frequency  as  a  variable  influencing  the 
resistivity. 

In  this  work  the  need  to  measure  potential  voltage  drops  radially 
across  the  sample  required  that  the  two-electrode  method  be  employed. 
Experiments  conducted  with  Berea  sandstone  samples  tested  with  the 
Simultaneous  Property  System  showed  that  the  problems  of  contact  resis¬ 
tance  and  frequency  dependence  could  be  minimized  by  coating  either  the 
electrode  or  sample  faces  with  conductive  silver  paint  and  by  selecting 
an  operating  frequency  above  700  hertz. 

Figure  Al-2  shows  calculated  resistivities  as  functions  of  fre¬ 
quency  obtained  using  the  two-electrode  method.  All  the  test  runs  show 
axial  resistivities  decreasing  with  increasing  frequency.  For  the  Berea 
samples  with  ends  coated  with  conductive  silver  paint,  resistivity  values 
become  essentially  independent  of  frequency  for  frequencies  greater  than 
about  700  hertz.  Similar  results  were  obtained  when  the  axial  heads  of 
the  SPS  were  silver  painted  instead  of  the  sample  ends. 

Resistivities  measured  using  a  two-electrode  circuit  were  checked 
for  accuracy  against  resistivities  obtained  from  a  four-electrode  con- 


114 


figuration.  Two-electrode  resistivity  values  were  measured  across  Berea 
samples  coated  on  the  ends  with  silver  paint.  Four-electrode  resistivity 
values  were  obtained  on  the  samples  using  rings  of  conductive  paint  on 
the  cylinder  circumference  as  measurement  probes.  Figure  Al-3  shows  the 
set-up. 

A  comparison  of  the  results  obtained  from  the  two-electrode  and 
four-electrode  methods  is  given  in  Figure  Al-4.  Resistivities  at  1,000 
hertz  calculated  from  potential  drops  measured  using  both  methods  are 
plotted  against  each  other.  The  relationship  is  linear  with  an  equation 
of  the  form  y  »  a  +  bx.  The  regression  coefficients  were  obtained  by 
least  squares.  Values  for  "a"  and  "b"  were  0.119  and  0.930  respectively. 
The  coefficient  of  correlation  for  the  regression  analysis  was  0.99. 

This  close  agreement  between  the  two-electrode  and  four-electrode  methods 
demonstrates  that  two-electrode  contact  resistance  can  be  minimized  and 
that  potential  drops  obtained  at  1,000  hertz  are  representative  of  the 
true  resistivity  of  the  brine  saturated  samples. 

In  order  to  apply  the  results  obtained  in  the  axial  resistivity 
measurements  to  radial  resistivity  experiments,  a  more  practical  and 
durable  electrode  than  silver  painted  stainless  steel  was  required.  Tests 
with  hard  industrial  chromium  plating  on  stainless  steel  radial  pass¬ 
through  shoes  were  successful  and,  as  with  the  axial  measurements,  gave 
resistivities  which  were  essentially  independent  of  frequency  when  fre¬ 
quencies  were  greater  than  700  hertz.  The  chromium  plated  passthrough 
shoes  were  molded  into  a  polyurethane  core  jacket,  and  the  resultant 
jacket-shoe  arrangement  was  used  to  contain  each  of  the  rock  samples. 


FIGURE  Al-2  FREQUENCY  DEPENDENCE  OF  RESISTIVITY  MEASUREMENTS  MADE  ON  TWO  BEREA  SANDSTONE  SAMPLES 
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4-  Electrode  Method 


FIGURE  Al-4  COMPARISON  OF  RESISTIVITIES  OBTAINED  FROM  TWO-ELECTRODE 
AND  FOUR-ELECTRODE  METHODS 
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B.  Application  of  Radial  Resistivity  Measurements 

For  axial  measurement  through  a  cylindrical  rock  sample  only 
the  current,  voltage  drop,  and  core  dimensions  are  necessary  for  re¬ 
sistivity  calculations.  However,  for  radial  measurement  through  a 
cylindrical  sample,  an  analog  to  the  permeability  formula  for  a  whole 
core  type  permeameter2 8  was  utilized.  Equation  3.9  is  such  an  analog: 
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The  geometric  and  decreased  current  flow  correction  factors  are 
based  on  ratio  values  of  passthrough  length  and  diameter  to  sample  length 
and  diameter.  For  this  work  Ke/Ke,  **  1.00  and  R^  “  0.94. 

Two  cut  samples  of  a  Berea  sandstone  rock  were  used  to  test  the 
reliability  of  resistivity  measurements  calculated  from  this  method. 
Berea,  if  properly  selected,  can  be  nearly  isotropic  and  homogeneous  for 
all  practical  purposes.  However,  for  comparison,  two  clean  Berea  cores 
were  cut,  one  perpendicular,  the  other  parallel  to  the  general  bedding 
plane  orientation.  Each  sample  was  saturated  with  brine,  inserted  into 
the  core  jacket,  placed  within  the  SPS,  and  subjected  to  current  flow 
axially  and  radially.  Figure  Al-5  depicts  the  flow  directions.  Axial 
resistivities  were  calculated  using 
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where  A  ■  area  perpendicular  to  current  flow  (cm2) 


Al.l 


Radial  resistivities  were  calculated  using  equation  3.9.  Results  showed 
that  axial  resistivities  of  Sample  A  were  within  2%  of  the  radial  re¬ 
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SAMPLE  A.  SAMPLE  B. 


FIGURE  Al-5  CUT  SAMPLES  OF  BEREA  SANDSTONE  USED  FOR  TESTING  RADIAL 
RESISTIVITY  MEASUREMENTS  (EQUATION  3.9) 
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slstlvltles  of  Sample  B.  No  trend  existed  and  resistivities  were  taken 
at  1,000  hertz.  The  results  verify  that  equation  3.9  provides  reliable 
resistivity  values. 

Rock  resistivity  data  in  its  basic  form  has  limited  use  because 
of  its  dependency  on  the  saturating  fluid.  The  formation  factor  concept 
developed  by  Archie11*  normalizes  the  rock  resistivity  and  renders  it  in¬ 
dependent  of  the  pore  fluid  resistivity.  The  formation  factor  as  defined 
by  equation  3.4  is  an  intrinsic  property  of  the  rock  at  a  specified 
effective  stress  state.  It  is  extremely  useful  so  far  as  it  can  be  re¬ 
lated  to  other  rock  parameters  such  as  porosity  and  permeability  as  dis¬ 
cussed  in  Chapter  III. 


C.  Pressure  Dependence  of  Resistivity  Measurements 

Pressure  effects  on  the  resistivity  of  formation  brines  usually 
are  not  significant.  Guyod16  explained  that  for  the  usual  formation 
waters  an  increase  in  pressure  from  0  to  5,000  psi  results  in  a  decrease 
in  resistivity  of  less  than  5 %  at  70°F. 

In  this  work,  however,  resistivity  measurements  were  taken  at 
various  effective  stress  states  and  at  pore  pressures  as  high  as  15,000 
psi.  For  this  reason  and  since  the  change  in  resistivity  of  the  rock 
with  effective  stress  was  the  variable  of  interest,  pressure  effects  on 
the  resistivity  of  the  saturating  fluid  had  to  be  taken  into  consideration. 

Quist  and  Marshall34  measured  electrical  conductances  of  aqueous 
sodium  chloride  solutions  from  0°  to  800SC  and  at  pressures  to  58,000  psi. 
However,  their  work  was  limited  to  solution  concentrations  up  to  only  0.? 
molar  (0.58%).  Therefore,  resistivities  as  functions  of  pressure  under 
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isothermal  conditions  for  NaCl  brines  of  6%  concentration,  which  was  the 
saturating  fluid  of  the  samples,  had  to  be  obtained  through  experiment. 

the  experimentation  consisted  of  placing  a  volume  of  6%  NaCl 
solution  in  a  Perspex,  thick-walled  cylinder  equipped  with  brass  end 
pieces  for  measuring  voltage  drops.  A  two-electrode  method  for  fluid 
resistivity  determination  was  employed.  Current  introduced  at  one  of 
the  brass  end  pieces  flowed  through  the  brine  and  out  the  other  end 
piece.  The  voltage  drop  across  the  brine  was  obtained  utilizing  the 
same  end  pieces  as  measurement  electrodes. 

Testing  occurred  at  various  pressures  up  to  15,000  psi.  To 
achieve  such  pressures,  the  Perspex  apparatus  was  placed  into  the  SPS. 

The  confining  pressure  generated  by  the  SPS  was  conveyed  to  the  brine 
solution  through  a  1/16  inch  tubing  piece  extended  with  a  12  inch  long, 
1/16  inch  ID  hose.  This  configuration  acted  merely  as  an  oil/water 
separator.  Figures  A1-6A  and  A1-6B  show  the  apparatus  between  the  heads 
of  the  SPS. 

Voltage  drops  across  the  brine  in  the  Perspex  cylinder  were 
measured  with  digital  volt  meters.  The  potential  drops  together  with  the 
current  flow  and  the  dimensions  of  the  cylinder  were  used  to  calculate 
the  brine  resistivity.  The  length  and  diameter  of  the  fluid  column  in 
the  Perspex  cylinder  changed  with  pressure  due  to  fluid  compressibility 
and  deformation  of  the  container.  The  dimensional  change  of  the  fluid 
column  was  corrected  for  in  the  resistivity  computations  by  measuring 
the  strain  of  the  hollow  cylinder  under  the  hydrostatic  stress  conditons. 
Rosette  strain  gauges  cemented  to  the  walls  of  the  cylinder  provided  the 
strain  measurements.  Since  Perspex  is  an  isotropic,  homogeneous  elastic 
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medium,  strain  can  be  expressed  as 


where  e  ■  strain  (inches /inch) 

v  ■  Poisson's  Ratio  (-  0.36  for  Perspex)35 
E  •  Young's  Modulus  (psi)  (■  4.5  x  105  psi  for  Perspex)35 
P  *  hydrostatic  pressure  (psi) 

Therefore,  the  internal  radius  and  length  of  the  brine  container  are 
given  by: 

r  -  ro(l  -  ^  P)  A1.3 

and 

L  -  L0(l  -  -y  2V  P)  A1.4 

where  ro  »  original  cylinder  internal  radius  (inches) 

Lo  *  original  cylinder  length  (inches) 

With  the  dimensions  of  the  fluid  column  corrected  for  deformation, 
the  resistivity  of  the  brine  at  various  pressures  was  calculated  using 
equation  Al.l.  The  results,  which  gave  resistivities  of  6?  NaCl  solution 
as  a  function  of  pressure  (temperature  constant) ,  are  shown  in  Figure 
Al-7.  Linearly  regressing  the  results  gave  the  following  relationship: 

Rg  -  0.1215  -  0.714P  A1.5 

where  Rg  ■  resistivity  of  the  brine  (ohm-meters) 

P  -  applied  pressure  (psi) 

Absolute  values  of  brine  resistivities  are  very  temperature  sen- 


sitive  as  shown  by  Guyod16.  It  was  found  during  this  experiment  that 
temperature  primarily  influenced  the  intercept  of  the  resistivity  verses 
pressure  line.  The  slope  of  the  line  remained  essentially  constant  for 
the  room  temperature  variations  which  occurred  during  the  rock  testing  ex¬ 
periments  . 

The  rock  resistivities  at  different  elevated  pore  pressures  were 
corrected  only  for  the  pressure  dependence  of  the  brine  resistivity.  To 
compensate  for  pressure,  formation  factors  were  calculated  using  brine 
resistivities  adjusted  for  pressure  effects. 

Temperature  was  not  a  factor  within  each  individual  rock  test 
since  each  test  was  conducted  under  essentially  isothermal  conditions. 

The  temperature  dependence  of  the  calculated  formation  factors  did 
arise,  however,  when  comparisons  were  made  between  individual  rock  tests. 
As  explained  above,  though,  the  slight  room  temperature  variations  merely 
changed  the  absolute  value  of  F  and  did  not  influence  relative  values. 
Therefore,  within  each  test  cycle,  variations  of  F  with  effective  stress 
could  be  compared  without  error. 


BRINE  RESISTIVITY,  Rg  (ohm  meters) 
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FIGURE  Al-7  GRAPH  OF  6Z  NaCl  SOLUTION  RESISTIVITY  AS  A  FUNCTION  OF 


PRESSURE 


APPENDIX  II 


Equipment  Calibration  and  Sensitivities 


Various  components  of  the  Simultaneous  Property  System  (SPS) 
required  calibration.  The  internal  axial  load  cell,  strain  gauges,  pore 
fluid  tubing,  and  wave  transducer  cell  heads  were  affected  by  the  pres¬ 
sure  and,  consequently,  were  calibrated  as  functions  of  the  confining, 
pore,  or  axial  pressure.  The  functions  were  essentially  linear  as  re¬ 
flected  in  the  calibration  equations  provided  in  Table  A2-1. 

The  waves  generated  and  received  by  the  transducer  cell  heads 
(Figure  3)  were  also  affected  by  the  epoxy  coatings  applied  to  the  ends 
of  the  sample  cores.  In  order  to  correct  measured  p-wave  propagation 
times  through  the  rock  samples  for  the  epoxy,  a  correction  factor  was 
determined  using  an  aluminum  cylinder.  The  cylinder  was  placed  in  the 
SPS,  and  p-wave  propagation  times  measured  at  various  pressures.  The 
same  measurements  were  then  taken  on  the  aluminum  cylinder  with  epoxy 
on  its  ends.  Comparisons  of  the  results  showed  the  p-wave  velocity 
through  the  epoxy  dependent  upon  the  confining  and  axial  pressures  applied. 
The  correction  factor  was  found  to  be  best  expressed  in  terms  of  the 
propagation  time  as  reflected  in  the  following  equation: 


where 


AT  -  (12.19  -  0.80  P.  -  0.64 

p  A  P,  t 


AT  ■  correction  factor  to  be  subtracted  from  p-wave 
p  propagation  times  measured  through  the  rock 
samples  (usee). 
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PA  *  axial  pressure  (kpsi) 

Pc  ■  confining  pressure  (kpsi) 

L„  -  total  thickness  of  the  epoxy  coatings  (Inches) 

The  pressures  and  sample  deformation  data  needed  for  the  cali¬ 
bration  equations  in  Table  A2-1  and  for  the  epoxy  correction  factor 
equation  (A2.1)  were  measured  with  digital  voltmeters.  The  recorded 
voltages  were  converted  to  useable  units  via  the  sensitivities  given  in 
Table  A2-2. 
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TABLE  A2-1 


Equipment  Calibrations  as  Functions  of  Pressure 


System  Designation 


pe  Device 


Calibration  Equation 


Axial  Pressure  Load  Cell  APA(kpsi)  *  (7.5  x  10  3)pc 

Sample  Deformation  Strain  Gauges  Ae(|-^)  *  (5.48  x  10  6)Pc 

Pore  Fluid  Volume3  Tubing  &  Valving3  AVw(cc)  *  (0.153)Pq 

p-Wave  Propagation  Transducer  Cell  Heads  AT  (ysec)  *  (5  x  10  3)P. 

p  A 


Notes:  1.)  All  Pressures  In  kpsi 

2. )  All  Calibration  Functions  Were  Subtracted  From  The 

Respective  Measured  Parameters 

3. )  The  Pore  Fluid  Volume  Change  With  Pressure  Was  A  Result 

Of  Not  Only  Tubing  And  Valving  Expansion  But  Also  Brine 
Compressibility31 ’ 36 
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TABLE  A2-2 


NOMENCLATURE 


A 

a 


C 

C' 


*  Area  perpendicular  to  current  flow 

*  Coefficient  in  the  generalized  Archie's 

*  Constant 
»  Constant 

■  Rock  bulk  compressibility 

■  Rock  matrix  compressibility 

*  Young's  Modulus 

*  Formation  factor 
=  Current  flow 

*  Geometric  correction  factor 

■  Permeability 
**  Axial  length 

»  Thickness  of  epoxy 

-  Original  Perspex  cylinder  length 

-  Cementation  factor 

*  Coefficients  of  the  cementation  factor 
«  Coefficients  of  the  formation  factor 

«  Pressure 

-  Axial  Pressure 

-  Confining  Pressure 

■  Effective  on  ess 

■  Pore  Pressure 
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equation 
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q 


“b 

Rf 


ro 


V 

V„ 


M 


eff 


v 

e 

£R 

G 

Z 

U 


=  Flow  rate 

=•  Resistivity  of  the  saturating  brine 

-  Correction  factor  for  decreased  flow  rate 

»  Resistivity  of  the  100%  brine  saturated  rock 
=  Internal  radius  of  Perspex  cylinder  under  pressure 
=  Original  Perspex  cylinder  internal  radius 

-  p-wave  propagation  time 
»  Voltage 

»  Bulk  volume 

«  Voltage  drop  across  resistor  measuring  contact  resistance 
=  Voltage  drop  of  resistivity  line  circuitry 

■  Voltage  drop  across  the  digital  voltmeter 
«  Initial  bulk  volume 

■  Pore  volume 

=  Voltage  drop  across  the  rock  sample 

■  Voltage  signal  from  signal  conditioner 
=  Water  volume 

•  Measured  p-wave  velocity  through  the  sample 

-  p-wave  velocity  through  the  matrix  material 
«  p-wave  velocity  through  the  brine 

-  Volumetric  strain  due  to  effective  stress 

-  Poisson's  Ratio 

-  Strain 

■  Radial  strain 

-  Axial  Strain 


Viscosity 


V, 
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$  »  Porosity 

<J> o  »  Initial  Porosity 

<pi  »  Porosity  by  water  volume  and  strain  measurement  method 

ip 2  ■  Porosity  by  water  volume  measurement  method 

<f>3  =  Porosity  by  volumetric  strain  measurement  method 

4> i»  »  Porosity  by  resistivity  measurement  method 

4>s  *  Porosity  by  Wyllie  Time  Average  method 

n  =  Coefficient  in  the  effective  stress  equation 

,s 
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